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Differential studies of elliptic flow are one of the most powerful tools in studying the initial conditions and 
dynamical evolution of heavy ion collisions. The comparison of data from Cu+Cu and Au+Au collisions taken 
with the PHOBOS experiment at RHIC provides new information on the interplay between initial geometry 
and initial particle density in determining the observed final state flow pattern. Studies from PHOBOS point 
to the importance of fluctuations in the initial state geometry for understanding the Cu+Cu data. We relate the 
elliptic flow data to the results of our model studies on initial state geometry fluctuations and discuss how we 
will perform measurements of event-by-event fluctuations in elliptic flow in Au+Au collisions. 
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I. ELLIPTIC FLOW 

Studies of the elliptic flow of produced particles via mea- 
surements of their azimuthal distribution are important probes 
of the dynamics of heavy ion collisions. In the collision of two 
nuclei with finite impact parameter, the almond-shaped over- 
lap region has an azimuthal spatial asymmetry. The asym- 
metrical shape of the source region can only be reflected in 
the azimuthal distribution of detected particles if the parti- 
cles significantly interact after their initial production. Thus, 
observation of azimuthal anisotropy in the outgoing particles 
is direct evidence of interactions between the produced par- 
ticles. In addition, these interactions must occur at relatively 
early times, since expansion of the source — even if uniform — 
would gradually erase the magnitude of the spatial asymmetry. 
Typically, hydrodynamical models are used to quantitatively 
relate a specific initial source shape and the distribution of 
emitted particles (2]. At top RHIC energy, ^/s^ = 200 GeV, 
the elliptic flow signal at midrapidity in Au+Au collisions is 
described under the assumption of a boost-invariant relativis- 
tic hydrodynamic fluid, indicating that there is early equilibra- 
tion in the colliding system J3]. 

The azimuthal anisotropy of the initial collision region can 
be characterized by the eccentricity (e) of the overlap region 
of the colliding nuclei in the transverse plane. The strength of 
the elliptic flow, vi, is commonly defined by the coefficient of 
the second harmonic in the Fourier expansion of the azimuthal 



particle distribution relative to the reaction plane, *Pr, such 
thatv 2 = (cos(24>-2'Fr)) 0. 

Recent hydrodynamical calculations J3l conclude that V2 
scales approximately with e for small e. Such calculations 
typically use smooth, event-averaged initial conditions, for 
which the initial azimuthal asymmetry is well described by 



the "standard" eccentricity, e s td 



where o 2 (a 2 ) is the 



variance of the participant nucleon distribution projected on 
the x (y) axis, taken to be along (perpendicular to) the impact 
parameter direction. 

PHOBOS has measured elliptic flow as a function of pseu- 
dorapidity, centrality, transverse momentum, center-of-mass 
energy J3-l7|] and, recently, nuclear species (H. In particu- 
lar, the measurements of elliptic flow as a function of central- 
ity provide information on how the azimuthal anisotropy of 
the initial collision region drives the azimuthal anisotropy in 
particle production. In fig. QJ, we show the centrality depen- 
dence of \>2 at mid-rapidity (|t|| < 1) for Cu+Cu and Au+Au 
at =62.4 and 200 GeV collision energies, as obtained 
from our hit-based analysis method iHH]. A substantial flow 
signal is measured in Cu+Cu at both energies, even for the 
most central events. This is quite surprising, as according to 
the initial anisotropy given by e st( j one expects that should 
approach zero as the collisions become more central, similar 
to what has been found for the Au+Au case. 
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FIG. 1: (a) v 2 (unsealed), (b) {e pa rt) and (c) v 2 /{£part) vs. N part 
for Cu+Cu and Au+Au collisions at ^/s m = 62.4 and 200 GeV. 1- 
a statistical errors (bars) and 90% C.L. systematic errors (bands) are 
shown. Data points and eccentricity calculation are from Refs. (aBl- 



II. ECCENTRICITY FLUCTUATIONS 

As a possible explanation for the large V2 signal in the small 
Cu+Cu system, we have argued that event-by-event fluctua- 
tions in the shape of the initial collision region may drive the 
elliptic flow iHlM]. For small systems or small transverse over- 
lap regions, fluctuations in the nucleon positions frequently 
create a situation where the minor axis of the overlap ellipse 
is not aligned with the impact parameter vector. These fluctu- 
ations are neglected in the definition of e st d- Instead, we have 
introduced the "participant eccentricity", 



(o?-a2)2 + 4o2 v 



^part 



(1) 



This definition accounts for the nucleon position fluctuations 
by quantifying the eccentricity event-by-event with respect to 
the minor axis of the overlap ellipse, in the frame, defined by 
"Ppart, that diagonalizes the ellipse. Note, a xy = (xy) — (x) (y), 
rj^ and a^, are the (co-)variances of the x and y participant nu- 
cleon position distributions expressed in the original frame, 
given by 'Pr. For a system with a large number of nucleons, 
the covariance term is comparatively small. Therefore, the av- 
erage values of e st( j and e part over many events are similar for 
all, but the most peripheral interactions for the Au+Au sys- 
tem. For the smaller Cu+Cu system, however, fluctuations in 



the nucleon positions become more important for all centrali- 
ties Jl. 

In fig. Hp, we show a Glauber model calculation of £ pa rt 
as a function of A^ pa rt for Cu+Cu and Au+Au. The collid- 
ing nuclei are built by randomly placing nucleons according 
to a Woods-Saxon distribution. Excluded volume effects are 
addressed by requiring a minimum inter-nucleon separation 
distance of 0.4 fm. As opposed to e s td, where averages are 
implicitly over participants and events, the variance expres- 
sions in e part are averaged event-by-event over participants, 
individually. To check how the event-by-event interpretation 
of the Glauber calculation depends on the external parameter 
settings, we varied a number of sources of systematic error, 
like the nuclear radius, nuclear skin depth, nucleon-nucleon 
inelastic cross-section Onn and minimum nucleon separation. 
Varying each specific parameter within reasonable limits, the 
individual contributions were added in quadrature to deter- 
mine the 90% confidence level errors shown in fig.[TJ). 

In order to compare the elliptic flow signals across nuclear 
species and with hydrodynamical predictions, it is important 
to scale out the difference in the initial asymmetry of the col- 
lision geometry. In fig. QJ, we show the eccentricity-scaled 
flow for Cu+Cu and Au+Au, vjj (e par t), as a function of cen- 
trality. The scaled data are very similar for both the Cu+Cu 
and Au+Au collision systems at the same number of partici- 
pants. It is important to note that the apparent scaling does not 
rely on the fine-tuning of the Glauber parameter settings, as is 
evident from the systematic errors, which are rather small (see 

fig.rip). 



III. ELLIPTIC FLOW FLUCTUATIONS 

If the proposed participant fluctuations are driving the value 
of V2 event-by-event, as suggested by our (v 2 ) analysis, then 
we expect them to contribute to event-by-event elliptic flow 
fluctuations. As mentioned above, ideal hydrodynamics leads 
to \'2 « e 01 . Assuming the same relation holds event-by- 
event, this would imply that a V2 /(v2) = <J e /(e), where o,, 2 
(rj £ ) is the standard deviation of the event-by-event distribu- 
tion of \'2 (e). Using our (v'2) data (fig. QJ) and the Monte 
Carlo Glauber simulation to obtain o £ (and (e)) for the par- 
ticipant eccentricity, we estimate a V2 to be about 2% for 
all centralities, except the most central Au+Au collisions at 
y/s NN = 200 GeV. This estimate leads to rather large relative 
fluctuations (a v , /(yi)) between 35 and 50%. It is important to 
note that these estimates neglect other sources of elliptic flow 
fluctuations. 

In the following, we will give a brief overview of the anal- 
ysis method that we have developed to measure flow fluctua- 
tions with PHOBOS in peripheral and semi-central Au+Au 
collisions. Details can be found in Ref. ifToll . We seek to 
discriminate known (statistical) from unknown (presumably 
dynamical) contributions to observed flow fluctuations. Ide- 

: a v 2 ,dyn + a v 2 ,staf 
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ally, they would add according to obs 
This relation holds if the average of the measurement, (Vj 
gives the true average in the data, (V2), and if the resolution 
of our method is independent of the true value. Neither of 
these conditions are fully satisfied in the event-by-event mea- 
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FIG. 2: Generated and reconstructed (v 2 ) (left) and <j,i 2 for AMPT 
vs. A'pan- The generated signal is extracted from a mixed-event anal- 
ysis based on the true Monte Carlo particle information. The recon- 
structed signal is averaged over results of running the complete flow 
fluctuation analysis in 10 vertex bins using the kernel created from 
the modified HIJING samples. The black lines represent the ±lo 
error on the mean, with the mean centered (not shown). 



surement of v 2 . Therefore, a more detailed knowledge of the 
response function is required. We define K(vf ,s ,v 2 ,n) as the 
distribution of the event-by-event observed elliptic flow, V2 bs , 
for events with constant input flow value, v 2 , and multiplicity, 
n. Assuming a set of events has an input v 2 distribution given 
by f(v 2 ), then g(v2 bs ), the distribution of v'2 bs , will be given 
according to 



K(vf s ,v 2 ,n)f(v 2 )N(n)dv 2 dn, 



(2) 



where N(n) is the multiplicity distribution of the events in the 
given set of events (centrality bin). Thus, our event-by-event 
elliptic flow fluctuation analysis consists of three steps: 

• Finding the observed v 2 distribution, g(i'2 bs ), for a set 
of events by an event-by-event measurement of v° bs . 

• Construction of the kernel, ^(v° bs , v 2 ,n), by studying 
the detector response for sets of constant (known) input 
value of v 2 and multiplicity n. 

• Calculating the true v 2 distribution, f(v 2 ), by finding a 
solution to eq. ||2). 

For the event-by-event measurement we use a maximum 
likelihood method, making use of all hit information from 
the multiplicity array to measure a single value, v£ , while 
allowing an efficient correction for the non-uniformities in 
the acceptance. We model the measured pseudorapidity de- 
pendence of v 2 according to V2(r|) = v 2 ■ (1 — |r||/6), with 



v'2 = v'2(0), which is known to describe our mean v 2 data rea- 
sonable well. We define the probability distribution function 
(PDF) of a particle to be emitted in the direction (t|,(|)) for an 
event with v 2 and reaction plane angle 0o as P(r|,0|v2,0o) = 
/?(r|)[l +2v2(t|)cos(2(|) — 20o)]- The normalization /?(r|) is 
constructed such that the PDF, folded with the PHOBOS ac- 
ceptance, yields the same value for different sets of parame- 
ters (i'2,(j)o). Maximizing nf=i-' 5 ( r li) < t ) il v 2j < t ) o) as a function of 
v 2 and 0o allows us to measure i'2 bs event-by-event. We de- 
termine the response function K(v° bs , v 2 ,n) in bins of v 2 and 
n using modified HIJING events. Flow of constant magnitude 
(v 2 ) with a flat reaction plane distribution (0o) is introduced 
into generated HIJING Au+Au events. This is achieved by re- 
distributing the resulting particles in each event in randomly 
according to 1 + 2v2(r|)cos(20 — 20o), using their generated 
T| positions. The modified HIJING events are run through 
GEANT to obtain the PHOBOS detector response. To finally 
extract f(v 2 ), we for now simply assume a Gaussian distribu- 
tion with iy 2 ) and a, 2 . For given values of (V2) and a, 2 , it is 
possible to take the integral in eq. (O to obtain the expected 
distribution, g^ v {v 2 *' s \ ( v 2),Gv 2 )- By comparing the expected 
and observed distributions, values for (V2) and a V2 are ob- 
tained from a minimized % 2 of the data to the expectation. As 
outlined in Ref. Oioll . the whole analysis procedure was veri- 
fied on similar HIJING events as used to construct the kernel, 
and found to successfully reconstruct the input fluctuations, 
provided (v 2 ) > 0.03. 

Here, we report on a study with fully simulated AMPT 
events in order to verify the analysis with a different type of 
"data" events than used to create the kernel, since in AMPT 
flow builds up dynamically. To minimize acceptance effects, 
the analysis is done in bins of collision vertex. Fig. [2] shows 
the averaged results obtained from the different vertex bins 
compared to the generated signal. Since the information about 
the generated v 2 in AMPT has not been readily available on an 
event-by-event basis, we extracted the generated (V2) and o,, 2 
from a mixed-event analysis based on the true Monte Carlo 
particle information. We conclude that the developed analysis 
chain is able to reconstruct the fluctuations to a satisfactory 
degree over a large range in centrality, in samples of "data" 
that are different from the ones used to construct the kernel. 
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